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Abstract 

We investigate fully nonlinear, non-Abelian excitations of quark-antiquark plasma, using relativistic 
fluid theory in cold plasma approximation. There are mainly three important nonlinearities, coming from 
various sources such as non-Abelian interactions of Yang-Mills (YM) fields, Wong's color dynamics and 
plasma nonlinearity, in our model. By neglecting nonlinearities due to plasma and color dynamics we get 
back the earlier results of Blaizot et. al. [I] . Similarly, by neglecting YM fields nonlinearity and plasma 
nonlinearity, it reduces to the model of Gupta et. al. [2], Thus we have the most general non-Abelian 
mode of quark-gluon plasma (QGP). Further, our model resembles the problem of propagation of laser 
beam through relativistic plasma [4] in the absence of all non-Abelian interactions. 
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1 Introduction 



QGP is a quasi-color-neutral gas of quarks and gluons which exhibits collective behaviour. It is expected 
to be formed in relativistic heavy ion collisions (RHICs) experiments, deep inside the neutron star and 
might have formed in early universe. A study of collective excitations of QGP are important to diagnose 
various parameters and signatures of QGP. It is also proposed that the chaotic collective modes of QGP 
gives an estimate of thermalization of QGP in RHICs [3]. From the extensive study of electrodynamics 
plasma, we know that there exists various linear and nonlinear excitations in it, governed by elcctrodynamic 
interactions which is an Abelian gauge theory. Here, in QGP, also we expect similar linear and nonlinear 
modes, but modified by the non-Abelian interactions, which itself is nonlinear. Therefore, in QGP, there 
is two types of nonlinear effects, one coming from usual plasma nonlinearity and other, from non-Abelian 
effects. Nonlinear solutions of non-Abelian or Yang-Mills (YM) theory is studied extensively by Matinyan 
et. al. without plasma, but with Higgs order phase [5]. Later, these studies were extended to QGP by 
Blaizot et. at, but without plasma nonlinearity and Wong's color dynamics [3] which is another non-Abelian 
effect in QGP. Various periodic, quasi-periodic, chaotic nonlinear modes and transition from order to chaos 
by plasma collective effects were studied. A study of stabilization of QCD vacuum instability by plasma 
collective modes were studied earlier in [7 . There is another group of work in these lines by Gupta et. at, 
but without the nonlinearity of YM fields and plasma nonlinearity. Nonlinear or Non-Abelian effects come 
from the Wong's color dynamics. Here, we present a fully nonlinear, non-Abelian excitations, including all 
nonlinearities: plasma nonlinearity, YM fields nonlinearity and color dynamics nonlinearity. 



the equation of motion, where m the mass, r the proper time, g the coupling constant, a the color index, u' 
4-velocity, which is also the fluid velocity in cold plasma limit, and G%" the field tensor, defined as, 



2 Relativistic fluid theory of QGP: 



The relativistic fluid set of equations, in cold plasma limit, is given by [5] 




(1) 



G?'=d»AZ-ff'AI£ + ge abc AZA v c 



(2) 
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in terms of 4- vector potentials A%. I a is the dynamical color charges which obey Wong's equation, 

^T- = -geabcU^Aublc ■ (3) 

CLT 

The vector potentials are obtained from the Yang-Mills field equation, 

duG? + ge abc A^ b G^ = J v a , (4) 

where J% is the 4-vector color current produced by various species in plasma with color charges, such as 
quarks, antiquarks and gluons. For simplicity, here in our analysis, we consider quark-antiquark plasma and 
the current density is given by 

r a = ]T nglaW, (5) 

species 

where n is the density of each species, determined by the continuity equation, 

d^nu^) = . (6) 

In general, these are very complicated, coupled, nonlinear equations to be solved and hence one goes for 
approximations, such as moving frame ansatz [3], space-homogeneous solutions, so on, to look for special 
solutions. Following Blaizot et. al. [T], let us consider the homogeneous solutions of our set of equations, 
few of them may be easily solved. The continuity equation for each species may be integrated and we get 

n(t)u°(t) — constant = iiqUq , (7) 

where no and Uq are the density and 0-component of fluid velocity at equilibrium. We also chose a gauge 
A° a = and the spatial part of the equation of motion may be easily integrated to get 



, _ gI a Aj 



(8) 



with the assumption that, at equilibrium, the plasma is at rest. The 0-component fluid velocity is given by 



u° = A /l + u?, (9) 



and hence Uq — 1. Similarly, the spatial part of the field equations gives 



# a + g 2 

and temperal component gives, 



{A{4)Ai - (A{Ai)Ai\ = g J2 noh ^ , (10) 



Y^ n oI a = eabcAlAi , (11) 



where dot means differentiation with respect to time. Finally, the color dynamics equation reduces to 

u i 

I a = ge ahc -^A\I c . (12) 

For further simplification, let us use hedgehog ansatz where the color directions are taken to be along the 
spatial direction and redefine variables as 

Y - 9!qA x i _ y _ gI A y2 z _ gI A z3 _ ^ 
TCI vci m ' ' 

and rescaling time and color charges as 

777 / 

t _> f and 4 -> f , (14) 
Jo Jo 

where Jo is introduced to normalize I a I a = 1, which is one of the constant of motion as can be seen from the 
equation for color dynamics. Further, from Eq. (jTTJ) , I a of second species (antiquarks) is opposite to that 
of first species (quarks) and hence I 2a = —I\a = —la- in terms of redefined variables, our simplified set of 
equations becomes, 

X + (Y 2 + Z 2 )X = -e ll , X , (15) 

y/1 + (I X X) 2 + {IyYf + (I Z Z) 2 ■ 



and 

r I y I z (Y 2 -Z 2 ) 



(16) 



y/l + (I x X)* + (I y Y)* + (I z Z)* 

and similar equations for y and z components which may be obtained by cyclic change among x, y and 

2oj 2 I 2 n 2 I 2 

z. The parameter e = J 2 " , where the plasma frequency to 2 = n ° 9 m . Above set of equations has two 
immediate constant of motions, namely, I a I a = 1 and 



(X 1 + Y z + Z z )/2 + {X Z Y Z + Y Z Z Z + Z z X z )/2 + e ^1 + (I X X) 2 + (I y Y) 2 + {I Z Z) 2 = E , (17) 

the energy. These approximate set of equations retains all important aspects of QGP such as YM nonlinearity, 
plasma nonlinearity and color dynamics nonlinearity. In the earlier calculations of Blaizot et. al. [T], the 
color dynamics is neglected and in [2], YM nonlinearity is dropped out. 

In order to extract the results of Blaizot et. al., let us assume that color charge I a are constant and then, 
Eq. (p~5|) reduces to, 

X + {Y 2 + Z 2 )X = -e\^ X , (18) 

3^1 + (X 2 + y 2 + Z 2 )/3 



where the square root term is the plasma nonlinearity, coming from the relativistic treatment just like in [?]. 
Further, expanding the plasma nonlinearity term up to 3 rd order in vector potential gives, 

X + (1--^)(Y 2 + Z 2 )X--^X 3 + ^X = , (19) 

which is similar to that of Blaizot et. al., except with few new terms containing a This new terms 

may lead to additional new features like chaotic scattering [S]. This model without these new additional 
term was studied extensively by Matinyan et. al. [5J and Blaizot et. al. pQ. 

Let us look for some other new solutions of our model Eq. (fT5|) . For example, a special solution with 
Z = leads to, 

X + Y 2 X = -e- ■ X (20) 
3 y/l + (X 2 + Y 2 )/3 

for X and a similar equations for Y with X and Y interchanged. It differs from similar work by Matinyan 
et. al., Blaizot et. al. because we have kept the plasma nonlinearity also. Our numerical study shows that 
the plasma nonlinearity enhances the chaos and therefore, increases the order to chaos transition parameter 
defined in [5 , which will be discussed later. 

Next let us look at another special solution with X = Y = Z of our general equation Eq. (fT5|) and we 

get 

o 1 X 

I + 2I 3 = -6- , (21) 

which describes a nonlinear oscillation. It is the more general nonlinear oscillation, including the plasma 
nonlinearity, than the elliptic functions Cn discussed in [5J[T]. It is easy to see that on neglecting the plasma 
nonlinearity, we get back Cn or Sn, depending on the strength of non-Abelian parameter compared to plasma 
frequency. It is interesting to note that above mode is an exact solution of QGP because for X = Y = Z, 
the color dynamics equation shows that color charges are constant. 

3 Results: 

The most general set of equations of our model are nonlinear, coupled equations and may not be easy to 
solve. So we have made an approximation, known as hedgehog ansatz and reduced the number of equations 
to be solved, but having all non-Abelian and nonlinear features. From these simplified set of equations, we 
may get the results of all other earlier works in this field. For example, in Fig. 1, we plotted the Poincare 
section of our model with the approximation that the dynamical color charges are constant and Z — 0, Eq. 



(HHJ). Figures 1(a) and 1(b) are for the system without plasma nonlinearity and found that the regular orbits 
seen in Fig. 1(a) for e = 5 disappears at the critical value of e = 2, Fig. 1(b) and hence chaotic. Similar 
figures with plasma nonlinearity shows changes from ordered orbit islands for e = 8.15, Fig. 1(c), into chaotic 
motion for e = 6. Therefore, the critical value of e for the order-to-chaos transition is higher with plasma 
nonlinearity. The chaos seen with e = 6.0 with plasma nonlinearity develops islands of ordered motion 
without plasma nonlinearity and we need smaller e (e = 2) to have chaos. Hence the plasma nonlinearity 
enhances the chaos which is an additional new feature compared to the results of Blaizot et. al. PQ . Another 
special solution of our model with I a — constant is X = Y = Z, Eq. (f2"Tj) . which is not an elliptic functions 
as in PP, but little more general nonlinear oscillation. 

Next, in Fig. 2, we plotted the general solutions of our model, as an example u x (Eq. (J5J) ) , with hedgehog 
ansatz for different values of e with the same initial conditions and we see that as the e decreases the system 
becomes more and more chaotic, which is, qualitatively, similar to the results of Gupta et. al.. For a large 
e, say, e = 100 (Fig.l (a)), the amplitude of oscillations is small and the YM nonlinearity and plasma 
nonlinearity may be negligible and hence it is just the Abelian oscillations, modulated by color dynamics. 
As e decreases, amplitude increases and all nonlincarities due to YM fields, color dynamics and plasma 
nonlinearity come into play and drive the system to chaotic motion as can be seen from Figs. 2(b) to 2(c) 
with intermittent oscillations. For e = 0, Fig. 2(d), the chaotic oscillations is mainly due to YM nonlinearity. 
Similar behaviour is also seen in the other components of velocity. 

4 Conclusions: 

We have studied fully nonlinear, non- Abelian excitations of quark- antiquark plasma using relativistic fluid 
theory. It exhibits new features like a special nonlinear oscillation, different from elliptic functions, and 
enhancement of chaos. Further, we have found that on neglecting color dynamics and plasma nonlinearity, 
we get back the results of Blaizot et. al. pQ and by neglecting YM fields nonlinearity and plasma nonlinearity, 
we obtain the results of Gupta et. al. [2J. Hence we have the most general nonlinear, non- Abelian modes 
of QGP. In general, all three nonlinearities are always there in the system. For small amplitude excitations 
(\X\, \Y\ and \Z\ are <C 1), YM field nonlinearity and plasma nonlinearity may be negligible like in [2J. 
At present, we don't know the appropriate limit to get the results of pQ, since the nonlinearity due color 
dynamics is always there irrespective of the strength of the excitations and the plasma nonlinearity also will 
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be there due to relativistic treatment. On neglecting all non-Abelian nonlinearities our model resembles the 
model of laser propagation through relativistic plasma [J] . 
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Figure 1: Poincare sections of our model (with Z = and I a = constant) without plasma nonlinearity (Figs, 
(a) for e = 5 and (b) for e = 2), and with plasma nonlinearity (Figs, (c) for e = 8.15 and (d) for e = 6). 



Figure 2: Exact numerical solutions of our model, as an example u x , for different values of e with the same 
initial conditions (Figs, (a) for e = 100, (b) for e = 20, (c) for e = 2 and (d) for e = 0). 



Figure Caption: 
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This figure "figl.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/0706.0771vl 



This figure "fig2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/0706.0771vl 



